function [Ydd,Ydq,Yqd,Yqq] = vectorize_mdl(im_dq_mdl,wscope)
%
[im_dq_mag, im_dq_phase] = bode(im_dq_mdl,wscope);
im_dq_phase = im_dq_phase/180*pi;

Ydd = reshape(im_dq_mag(1,1,:),[],1).*...
    exp(1i*reshape(im_dq_phase(1,1,:),[],1));

Ydq = reshape(im_dq_mag(1,2,:),[],1).*...
    exp(1i*reshape(im_dq_phase(1,2,:),[],1));

Yqd = reshape(im_dq_mag(2,1,:),[],1).*...
    exp(1i*reshape(im_dq_phase(2,1,:),[],1));

Yqq = reshape(im_dq_mag(2,2,:),[],1).*...
    exp(1i*reshape(im_dq_phase(2,2,:),[],1));
